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ABSTRACT 

We use a new contour-based map analysis technique to measure the mass and size of molecular 
cloud fragments continuously over a wide range of spatial scales (0.05 < r/pc < 10), i.e., from the 
scale of dense cores to those of entire clouds. The present paper presents the method via a detailed 
exploration of the Perseus Molecular Cloud. Dust extinction and emission data are combined to yield 
reliable scale-dependent measurements of mass. 

This scale-independent analysis approach is useful for several reasons. First, it provides a more 
comprehensive characterization of a map (i.e., not biased towards a particular spatial scale). Such a 
lack of bias is extremely useful for the joint analysis of many data sets taken with different spatial 
resolution. This includes comparisons between different cloud complexes. Second, the multi-scale 
mass-size data constitutes a unique resource to derive slopes of mass-size laws (via power-law fits). 
Such slopes provide singular constraints on large-scale density gradients in clouds. 
Subject headings: ISM: clouds; methods: data analysis; stars: formation 



1. INTRODUCTION 

Some of the most fundamental properties of molec- 
ular clouds are the mass and size of these clouds and 
their substructure. Today, these properties are well con- 
strained: we know the masses and siz es of dense cores in 
molecular clouds (5, 0-1 PC size; e.g. iMotte et al.lll998l 
Johnstone et al.ll2000L lHatchell et al J I2005L lEnoch et al.1 
20071 ). and those of clumps (som e 0.1 pc) and clouds 



(^, 10 Pc) containing the co res (e.g.. I Williams et al.lll99l 
lCambresvll999llKirk et al.l2006t see lWilliams et al.l2000l 

for definitions of cores, clumps, and clouds). We do, 
however, not know much about the relation between the 
masses and sizes of cores, clumps, and clouds: tradition- 
ally, every domain is characterized and analyzed sepa- 
rately. As a result, it is still not known how the core 
densities (and thus star-formation properties) relate to 
the state of the surrounding cloud. 

In principle, the relation between the mass in cloud 
structure at large an d small s patia l scales is described by 
mass-size relations. iLarsonl (|1981f ) presented one of the 
first studies of such relations. He concluded (in his Eq. 
5) that the mass contained within the radius r obeys a 
power-law, 



m(r) = 460 M (r/pc) 



1.9 



(1) 



Most subsequent work refers to this relation as "Lar- 
son's 3 rd law", and replaces the origina l result with 
m(r) oc r 2 (e.g.. IMcKee fc Ostrikerl 12007?) . This "law 
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of constant column density" (with respect to scale, 
r) is now considered one of the fundamental prop- 
erties of molecular cloud structure (e.g., reviews b; 



[ Ballesteros -Parcde s etjffll2007l IMcKee fc Ostrikerll200 
iBergin fc Tafallal 120071 ). This relation has, however, 
never been re-examined comprehensively on the basis 
of up-to-date data. It is, e.g., not clear whether re- 
cent dust extinction and emission work is consistent with 
m(r) oc r 2 . 

Further, the limitations of available s tructure identi- 
ficati on schemes (such as CLUMPFIND; Williams et al. 
Il994f ) forced past studies to break cloud structure maps 
up into discrete fragments. These fragments typically 
have a size slightly larger than the map resolution. As 
a consequence, the cloud structure is only probed in 
a narrow spatial domain; the largest spatial features 
in a given map are, for example, usually not charac- 
terized. Today, approaches permitting automatic ex- 
amination of_a_TOntinuo_us range o f spatial scales are 
available. Rosolowsky et al. (2008b), in particular, pro- 
vide software for such studies (their dendrogram anal- 
ysi s); our work woul d be impossible without the work 
by iRos olowsk v et ail Such software permits derivation 
of spatially more comprehensive mass-size relations than 
possible in the past. 

In this series of papers, we combine contemporary col- 
umn density observations of high sensitivity with a new 
data analysis technique to examine the mass-size rela- 
tion in molecular clouds for a continuous range of spatial 
scales of order 0.01 to 10 p c. We rely on extin ction maps 
of molecular clouds (h ere: iRidge et all 1200 6). as well as 
maps of dust e mission (lEnoch et al.ll2006l) . Following the 
terminology of iPeretto fc Fuller! (|2009f ) . we define "cloud 
fragments" in the maps as regions enclosed by a contin- 
uous column density contour, and derive their mass and 
size at various contour le vels. This is implemented u sing 
algorithms introduced by IRosolows kv et al.l (|2008bf ). 

The first two papers in this series establish our analysis 
approach (part I, the present paper) and explore several 



2 



Kauffmann et al. 



clouds in the solar neighborhood (< 500 pc; part II). 
Section [2] of the present paper describes our cloud frag- 
ment extraction and characterization scheme. In Sec. [3] 
we provide a first idea how basic physical properties and 
observational limitations affect the mass-size measure- 
ments. This includes a comparison to results obtained 
using the CLUMPFIND algorithm. A detailed discussion 
of analysis uncertainties is presented in Sec.[4j These are 
explored using data for the Perseus Molecular Cloud. We 
also explain how dust emission and extinction data can 
be combined for a given cloud. 

Section [5.11 briefly describes how the new map analysis 
scheme might help to advance star formation research. 
As we describe there, it will help to jointly analyse data 
taken at different spatial resolution. This is a key fea- 
ture in the age of multi- wavelength and multi-resolution 
studies. We conclude with a summary in Sec. 15.21 

2. METHOD 

2.1. Processing of Contour Maps 

Consider some column density map containing a num- 
ber of local maxima, as sketched in Fig. [T] (a). In 
this map, cloud fragments can be identified as struc- 
tures bound by some continuous contour. To charac- 
terize these, we pick one of the column density maxima, 
and measure the mass, m, and area, A, contained within 
each column density contour containing this peak. Here, 
we use the effective radius, 

r = (A/n) 1 /\ (2) 

to quantify A. By following the trends from contour to 
contour, as shown in Fig.[T](b), one can construct a mass- 
size relation for every cloud fragment. Strictly speaking, 
we do thus construct mass-area diagrams. In practice, 
however, r is arguably a more intuitive variable than A. 

Some contours may contain several local maxima. 
They represent composite fragments. To give an example 
based on Fig.[TJ the region bound by the dark-green con- 
tour consists of the two regions marked in lighter shades 
of green, and the magenta contour contains the fragments 
marked in yellow and dark-green. Merging of two frag- 
ments occurs at the first column density contour contain- 
ing both of the fragments. In our analysis, we enforce 
that only two fragments can merge at a time. 

During a merger, mass and size jump discontinuously 
from the pre-merger situation, m, and fj, to their post- 
merger value, m — mx + tob and r = (r\ + rf,) 1 / 2 . 
This yields gaps in the mass-size relation (Fig. ft] [b]). 
The merger information is preserved during processing, 
so that it is possible to look up which fragments are con- 
tained within others. 

Our map characterization s cheme is thus clo s ely re - 
lated to the one employed by iPeretto k Fulled (|2009f) . 
Like us, these authors measure sizes and masses for 
regions bound by lines of constant column density. 
The schemes differ in the number of contours consid- 
ered: we use a v ery large number (10 2 to 10 3 ), where 
IPeretto k Fulled only consider two contours per cloud. 

In practice, we us e the dendrogram ( i.e., tre e analysis) 
code presented by Ro solowskv et al.l ((20 08b) for auto- 
matic processing of the maps. A minimum significant 
contour has to be set for every region; emission below 



this limit is not characterized here. It is further neces- 
sary to specify a minimum contrast between peaks, in 
order to identify significant central maxima for the ob- 
jects. These threshold column densities and contrasts are 
listed separately for every region in the following. The 
identified maxima are required to be spaced by more than 
one spatial resolution element. Again, this parameter is 
noted separately for every map. In this work, all sources 
with an effective diameter (i.e., 2r) smaller than twice 
the map resolution are rejected and are treated as parts 
of enveloping objects. 

In the terminology of lRosolowskv et all our mass mea- 
surement approach (i.e., integrate column density within 
some contour) corresponds to their "bijection paradigm" . 
Such mass measurements, e.g. made towards a dense 
core, always give a sum over several spatially overlap- 
ping components (e.g., some fraction of the dense core, 
plus some fraction of its envelope). This is not a prob- 
lem, if properly taken into account in the later analysis. 
Ot her choices are pos sible (e.g., the "clipping paradigm" 
of iRosolowskv et all ) , but they are less intuitive and 
even harder to model. Also, most previously existing 
data has been published effectively adopting the bijec- 
tion paradigm. 

2.2. Mass Estimates 

As a first example, below we present a mass-size 
analysis of Perseus, based on column densities derived 
from 2MASS-deriy ed extinction data. As explained by 
iRidge et al.l ((2006), the map is derived in terms of mag- 
nitudes of visual extinction, Ay- We convert this to H2 
column densities using the relation 

N H2 = 9.4 x 10 20 cirT 2 (A v /mag) (3) 

(Bohl in et al.l lT978). Mass surface densities, E, can then 
be derived as 

E = ^h 2 m H Nn 2 , (4) 

where /iu 2 = 2.8 i s the mean mol e cular weight 
per H2 molecule ( Kauffm ann et all I2008T ) and 
mjj is the weight of the hydrogen molecule. In 
practice, E = 0.047 gem" 2 (A^/IO 22 cm~ 2 ) = 
226 M Q pc- 2 (A^h 2 /10 22 cm^ 2 ), or A v = 
227.5 mag (E/gcm -2 ). The mass is then derived 
by integrating the mass surface density, m — f E dA. 
For t his we adopt a Perseus distance of 260 pc (|Cernisl 

[1993T1 

iGoodman et all (|2009bl ) present a comparison of col- 
umn density tracers for the Perseus region. After study- 
ing column density maps based on dust extinction (from 
2MASS data), dust emission (from IRAS imaging), and 
line emission data (from large field 13 CO [1-0] mapping), 
they conclude that extinction-based maps provide the 
best available information on a cloud's spatial mass dis- 
tribution. Extinction-based column densities deviate by 
~ 25% from those derived from other tracers (after re- 
moving global offsets between estimates, e.g. due to the 
choice of dust opacities). The true column density is sup- 
posedly in between these estimates. If every tracer has a 
similar scatter with respect to the true column density, 
this scatter would then be ~ 25%/2*/ 2 w 18% for all 
tracers. Extinction-based estimates of the column den- 
sity do probably deviate by a lower amount from the true 
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Effective Radius [pc] 



Fig. 1. — Fundamental concept of mass-size measurements. Panel (a) shows an example column density map for the Perseus molecular 
cloud, where labels refer to individual star-forming regions. Such maps usually contain several local maxima (numbered crosses in map). 
We pick one of these maxima, and draw contours at constant column density around this peak (map inset). For each contour, we measure 
mass and size. These measurements can then be placed in a mass-size plot, as marked by crosses in panel (b). When a map contains 
several maxima, it can be divided into cloud fragments where contours just contain a single maximum (light green and yellow boundaries 
in panel [a]), two of these (dark green boundary), or even more (magenta contour). In this sense, mass-size measurements for contours 
contained in one of these areas are related. These relations are highlighted by colored lines in panel (b); the color refers to the map region 
from where the measurements are taken, and numbers indicate the central column density peak of each fragment. When two fragments 
blend into a combined one that contains both, mass and size measurements jump discontinuously from those of the individual fragments 
to those for the combined one. These jumps are indicated by dotted lines. 



value. Here, we thus adopt a systematic uncertainty of 
< 15% for extinction-based mass estimates in Perseus. 

Young stars embedded in the clouds can further bias 
extinction observations, given their red intrinsic colors. 
This bias is particularly significant towards clusters, such 
as NGC1333 and IC348 in Perseus. We do not exclude 
these regions from our study, but measurements towards 
the clusters should be interpreted with particular cau- 
tion. 

2.3. Example Map 

Figure [5] shows the aforementioned extinction map for 
Perseus. The two most prominent star-forming regions in 
Perseus, NGC1333 and IC348, manifest as extended ex- 
tinction structures in this map. Fragments of very small 
size (e.g., < 0.1 pc), like dense cores around individual 
young stellar objects, are not visible in the map, due 
to its too poor spatial resolution (5', corresponding to 
w 0.4 pc). Section T4.2I shows how such structures can 
still be bootstrapped into mass-size studies. 

Figure [3] presents the mass-size data derived from the 
Perseus extinction map. This diagram reveals that the 
Perseus complex has a total (effective) radius of 8 pc and 
a total mass of 1 x 10 4 M . The properties from con- 
tours containing the two major stellar clusters, NGC1333 
and IC348, are highlighted by bold black lines. As one 
may naively expect, the cloud fragments enclosing these 
clusters are, at given radius, the most massive fragments 
within the cloud complex. This analysis also reveals frag- 
ments that are, again at given radius, much less massive 




Fig. 2. — Example column density map for Perseus, presented in 
terms of visual extinction, Ay. The map is taken from Rid ge et al.l 
(2006). Contours are drawn in steps of 1 mag, starting at 2 mag. 
Labels indicate the rough position of the star forming regions 
NGC1333 and IC348. 



than the regions containing the clusters. As discussed 
in the next paragraphs, the column density sensitivity 
of the map sets a radius-dependent lower limit to the 
masses that can be detected in a given map. Fragments 
of a mass much lower than those shown here may thus 
well exist in Perseus. Eventually, at sufficiently large ra- 
dius, most of the contour-bound objects found merge into 
a single fragment containing essentially all of the Perseus 
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Fig. 3. — Mass-size relation for the Perseus molecular cloud. 
Compared to Fig. [T] we only keep the lines connecting related 
measurements. The top panel describes the data in detail. Bold 
solid lines indicate data for fragments either containing NGC1333 
or IC348. Light solid lines show observations for regions not as- 
sociated with these clusters. The dashed line gives the sensitivity 
limit of the analysis. In the bottom panel we present the data in 
the context of reference mass-size relations. Dotted lines indicate 
mean H2 column densities, (JVh 2 ), respectively the mass-size re- 
lation for a singular hydrostatic equilibrium sphere suppo rted by 
isothermal pressure from gas at 10 K temperature (see Sec. I3.ll for 
both). 

molecular cloud. 

In the processing of the map, we have used a minimum 
threshold extinction of 2 mag. T his limit is much l arger 
than the noise level of 0.4 mag (|Ridge et al.l 120061 ) and 



rejects the signal from background extinction structures 
not related to the cloud. Correspondingly, we cannot 
detect fragments with a mean extinction below 2 mag. 
Given the relation between mass, size and mean column 
density (Eq. [5]), this minimum column density sets a 
lower limit to the detectable mass. Substitution of the 
minimum column density, -ZVj^min (in Eq. [SJ note that 
(-^h 2 ) > Ah 2 min), gives sensitivity limits like the one 
shown in Fig. j§] The minimum contrast between peaks 
is required to be the noise level times a factor 3, i.e. 
1.2 mag. The spatial resolution of the map is 5'; we do 
therefore require the maxima to be separated by at least 
3 pixels (7.'5), and objects with a radius smaller than 5' 
are rejected as being unphysical. 

3. PROPERTIES OF MASS-SIZE DATA 

Some mathematical and physical laws governing the 
properties of mass-size data have to be heeded, if a mean- 
ingful interpretation of the observations is desired. Here 
we list the most fundamental of these. 

3.1. Reference Relations 

A number of reference mass-size relations can help to 
navigate within the observational data more intuitively. 
They are in part derived assuming a spherical geome- 
try for cloud fragments. The assumption of spherical 
symmetry may not be appropriate, though. This caveat 
should be kept in mind when using the following refer- 
ence relations. 

Mass and size measurements can be used to calculate 
the mean mass surface (or column) density of a frag- 
ment, (£) = m(r)/A(r). Conversely, one can draw lines 
of constant mass surface density, m(r) = (£) nr 2 . Con- 
version to column density, and substitution of the afore- 
mentioned constants, yields 

m(r) = 71 Mq «A^h 2 >/[10 21 cm" 2 ]) (r/pc) 2 . (5) 

These lines are drawn in most mass-size plots presented 
here (e.g., Fig. Equation (|5|) implies one of the most 
important properties of mass-size data: since column 
density decreases with increasing radiu£0, the observed 
mass-size relations must be flatter than m(r) oc r 2 . 

Equation ([5]) can actually be used to transf orm our re- 
sults into the analytical diagram presented bv lTanl (|2007t 
column- densi t y vs. mass). Our study goes beyond the 
work by [Tan (2007) in that it systematically populates 
the parameter space with observational data. 

For spherical clouds, the mean density is (g) = 
ro(r)/(4/3 7rr 3 ). The corresponding mass-size relation 
reads m(r) — 4/3 ir (g) r 3 , or 

m(r) - 282 M Q «n H2 )/[100 cm" 3 ]) (r/pc) 3 , (6) 

where we substitute the density of H2 molecules, uh 2 — 
g/(p-H2 m a)- 

In part II of this series, we shall study spherical power- 
law density profiles, g(s) cx s~ k (where s is the radius), 
as models for the observed mass-size relations. We show 
that 

g(s) oc s~ k m(r) oc r 3 ~ fe . (7) 

4 By definition, this is always the case for our source character- 
ization scheme. We start off from a column density peak and then 
consider, by design, regions that increase with size when lowering 
the threshold column density. 
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The slope of the mass-size relation is therefore related 
to the slope of the density law. A density profile 
p(s) oc s~ 2 is often adopted to describe dense cores (see 
iDapp fe Basull2009l for a discussion). This gives a mass- 
size relation m{r) oc r. In hydrostatic spheres supported 
by isothermal pressure, mass, size, and gas temperature 
are related by 

(see iKauffmann et al.1 120081 Eq. 13). For gas tempera- 
tures T g = 10 K, one obtains the model relation drawn 
in most mass-size diagrams of this paper (e.g., Fig. [3]). 

We stress that we obtain two-dimensional mass-size 
relations from column density maps. These are related 
to, but not identical with, mass-size laws obtained from 
three-dimensional density map s. This is illustrat ed by 
the experiments conducted by iShettv et al.l ([2010, sub- 
mitted) who use the fragment identification technique 
also used by us. Their analysis is based on three- 
dimensional numerical simulations of turbulent clouds. 
As part of their experiments, they fit power-laws (sim- 
ilar to Eq. [T]) to their mass-size data. For their par- 
ticular set of simulations, the exponent derived in this 
fashion is similar to the number of dimensions used for 
mass measurements (i.e., 3 when based on density, and 
2 when based on column density). This underlines that 
the number of dimensions considered has to be kept in 
mind. Note, though, that these details do not compro- 
mise mass-size measurements as a tool for cloud structure 
analysis. Observed mass-size laws unambiguously sum- 
marize actual cloud structure. Only their interpretation 
is sensitive to the assumed geometry. 

3.2. Relation to CLUMP FIND -like Results 

The approach chosen here constitutes one of several 
possible choices to measure the mass and size of ob- 
jects in maps. Another p opular approach is to use the 
CLUMPFIND algorithm (jWilliams et al.|[l99l . This 
method uses contours to break emission up into sev- 
eral objects, just as done here. CLUMPFIND-extracted 
boundaries do, however, not necessarily follow contours. 
This is a major difference to our method, where ob- 
jects are always bound by some column density contour. 
Based on this fact alone, CLUMPFIND a nd our approach 
will thu s extract very different objects. iGoodman et al.1 
(2009a) illustrate this problem comprehensively. 

Further, CLUMPFIND does not allow for hierarchical 
structure, i.e., fragments containing other fragments. For 
example, CLUMPFIND will never determine the integral 
properties of the entire cloud, since the cloud is usually 
broken up into many independent fragments. The rela- 
tion o f CLUMPFIND-resu lts to cloud hierarchy is stud- 
ied by EiidaiFil] (gOOl) . 

Figure [4] compares mass and size measurements from 
CLUMPFIND to those from our approach^- Both char- 
acterizations are based on our Perseus extinction map. 
To initiate CLUMPFIND, we choose contours as close as 
possible to the para meters used in our contour-based seg- 
mentation ( Sec. 12. 3]) ; the contour spacing is 1.2 mag, with 
the lowest contour at 2 mag. Still, this CLUMPFIND 

5 The CLUMPFIND results were kindly provided by J.E. Pineda. 
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Fig. 4. — CLUMPFIND (crosses) and dendrogram mass-size re- 
sults (lines) for the extinction map presented in Fig. [2] The two 
approaches yield broadly similar, but not directly related results. 

segmentation yields objects that are only remotely re- 
lated to the regions extracted by our method. To start 
with, CLUMPFIND extracts 52 peaks enclosed by in- 
dividual clumps, where our method only identifies 17 
significant peaks. This discrepancy is a consequence of 
CLUMPFIND 's relaxed peak identification criterion: ev- 
ery local peak encircled by a continuous contour is con- 
sidered significant, independent of the depth of the col- 
umn density dip towards the next peak. There is thus 
no well-defined correspondence between regions identi- 
fied by the different algorithms. 

To characterize the relation between objects extracted 
by different methods in some sense, we inspect the struc- 
ture around the 17 significant column density peaks 
found by our method. For every such peak, we iden- 
tify the CLUMPFIND object containing this peak. The 
mass of the latter clump can then be compared to our 
results, taken at the clump's radius. Because of merging 
of objects, our method does not provide a mass mea- 
surement for every possible radius (consider, e.g., the 
evolution of mass-size measurements starting from peak 
3 in Fig.QJb]). For those cases where mass measurements 
exist for the CLUMPFIND-derived radius, we find that 
CLUMPFIND gives masses of order 55% to 95% of the 
masses derived by our approach. 

In summary, mass and size measurements from 
CLUMPFIND are thus broadly compatible with our re- 
sults. By this we mean that the CLUMPFIND-derived 
mass-size measurements reside in the space spanned by 
our own measurements. There is, however, no good cor- 
respondence on a detailed level. 

4. PERSEUS IN DETAIL 
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Given the tools derived above, it is now possible to 
characterize the mass-size relation for Perseus. While do- 
ing so, we also evaluate uncertainties affecting our anal- 
ysis. In a hnal step, we extend the analysis to clouds 
other than Perseus. 

4.1. Mass Uncertainty 

To explore the impact of noise on mass measurements, 
we run Monte-Carlo experiments in which we add arti- 
ficial Gaussian noise of root-mean-square (RMS) ampli- 
tude a(mi) to the map before structure characterization 
(where is the mass per pixel). Such trials are pre- 
sented in Fig. EJa). Comparison of the derived mass-size 
data to the one from the original map does then reveal 
the impact of noise. In our experiments, we find that 

m m rbeam 

is an upper limit to the noise-induced mass changes, 
a(m). In this, rbeam is the beam radius and the nu- 
merical constant does slightly depend on the number of 
pixels per beam. 

Equation ([9]) follows from Gaussian error propagation 
of m — Sjrrij. The numerical constant is increased by 
a factor 3, though, to provide a strict upper limit to 
uncertainties. To test for non-Gaussian sources of error, 
we validate Eq. (J9j> by comparing the original data with 
those derived from maps with additional noise. For any 
given structural branch with r < 2 pc, Eq. ([§]) is indeed 
found to set an upper limit to the mass deviation at given 
radius, following the noise experiments depicted in Fig. 
[5ja). For larger radii, a small (but significant) additional 
uncertainty of ~ 1% has to be included to capture non- 
Gaussian sources of error. 

In our study, the relative uncertainty of column density 
estimates (i.e., a(vni)jmi) is < 1/5 in characterized re- 
gions (given extinction noise levels and selection thresh- 
olds of 0.4 mag and 2.0 mag, respectively). Further, 
the mass is at least as large as the one at the sensi- 
tivity limit (Fig. [Ha]), m > 7rm !iinin (r/r bcam ) 2 (in our 
map, the beam contains, by area, 7r pixels). Substitu- 
tion of these parameters in Eq. © give uncertainties 
of 19% for the smallest extracted features (for which 
r = 2rbeam) located just at the sensitivity limit. Since 
a(m)/m cx (-/Vh 2 ) -1 r -1 (via substitution of Eq. [5]), 
larger regions well above the detection threshold will suf- 
fer lesser uncertainties, <C 10%. 

In later papers, we will compare the properties of 
clouds that are located at different distances. Since the 
angular resolution of the observations is about constant, 
an increase in distance implies a decrease in physical reso- 
lution. To first order, this decrease in physical resolution 
corresponds to smoothing of the map. We explore the 
impact of smoothing by extracting mass-size data from 
smoothed maps of the Perseus cloud, as demonstrated in 
Fig. Eft). _ 

Smoothing does effectively mean that mass is trans- 
ferred to larger scales. For fixed radius, smoothing does 
thus imply a reduction in mass. To estimate this re- 
duction, consider a region of effective radius r. After 
smoothing, the mass contained in this region will be 
smeared our over an area of radius r' ~ r- m i + r\. In 
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Fig. 5. — Impact of noise (top) and limited spatial resolution 
(bottom) on mass-size measurements. Note that the size and mass 
scale differ between the panels. Noise and resolution are explored 
using artificial maps created by adding noise to (top), respec tively 
smoothing of (bottom), the observed data for Perseus (Sec. I4.U . 
Differences between the input data and the properties of the artifi- 
cial maps are due to these biases. The original input data is drawn 
in black, yellow lines in the top panel show mass-size relations from 
maps with additional noise similar to the observed one (0.4 mag). 
Red and yellow lines in the bottom panel present mass-size rela- 
tions for maps with a resolution worsened to rb ca m,f/f*bcam,i = 3/2 
and 5/2, respectively (and features with r < 2rb ea mf removed). 
At gi ven radius, both biases just induce moderate changes in mass 
( Sec. ED , 

this, r k = (rt am)f - rg oam , i ) 1/2 /(ln(2)) 1 / 2 is the effec- 
tive radius of the smoothing kernel required to go from 
the initial to the final beam radius of the map during 
smoothing, r bcam ,i -)• r hcam j. Very approximately, after 
smoothing the mass retained in the initial area will be of 
order of the fraction of the initial area to the one after 
smoothing, r 2 / (r + r^) 2 . 
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After rearrangement (and using <r), we find that 
the smoothing-induced mass bias should obey a relation 
of the form 



Am 



0.2 to 0.3 ■ 



1 



1 



(1 + r k /r) 5 



(10) 



To allow for more realistic source geometries, the nu- 
merical constant must be derived from experiments with 
actual data. We use the smoothing experiments shown 
in Fig. EJb) for this purpose. With these parameters, 
Eq. (fTU|) describes the relative mass bias for all structure 
branches in our example map. 

In our study, we only consider regions with radii larger 
than twice the effective beam radius, r > n^am- Sub- 
stitution of this (with rbeam,i = 0, since the true dis- 
tribution on the sky has infinite resolution) in Eq. (jlOl) 
shows that the mass of regions with a radius similar to 
the beam diameter may be underestimated by not more 
that 20%, and by less for larger fragments. 

4.2. Combining Dust Emission and Extinction Data 

Our Perseus extinction map has a spatial resolution 
of 5' (0.4 pc). It cannot be used to characterize much 
smaller fragments, such as dense cores of < 0.1 pc size. 
This limitation can, however, be overcome when includ- 
ing mass-size data from additional data sets. 

4.2.1. Basic Concept 

For Perseus, lEnoch et al.l ()2006l ) present a Bolocam 
map of dust emission in Perseus at 1.1 mm wave length. 
The beam width at half sensitivity is 31", corresponding 
to 0.04 pc. At these wave lengths, the continuum emis- 
sion of dust at temperature is optically thin, and so 
the observed dust continuum emission intensity, I v , is a 
direct measure of the column density along a given line 
of sight, 

I v =MH 2 TOH^iVH 2 -B ly (Td) (11) 

(see iKauffmann et al.l 120081 for full details), where B v is 
the Planck function, and the dust opacity, K v , IS eval- 
uated per total gas mass. Column density maps from 
dust emission can then be contoured and analyzed as de- 
scribed before (Sec. [2]) to derive mass-size data for the 
fragments in the dust emission map. 

Figure [5] combines the Bolocam-derived mass-size mea- 
surements with the extinction-based ones. Spatially, 
they are separated by a gap, since Bolocam (< 2') and 
2MASS extinction maps (> 5') probe different spatial 
scales. The extinction-identified fragments do, however, 
contain the Bolocam-derived ones. For example, shad- 
ing in Fig.[B]highlights Bolocam-detected fragments con- 
tained in the extinction peak harboring NGC1333. 

4.2.2. Relative Calibration of Dust Opacities 

iBianchi et all ((20031 ) studied B68 to examine the re- 
lation of mass estimates from dust extinction and emis- 
sion. They find relative difference s by factors of 1 to ~ 
2.35, w hen adopting opacities for lOssenkopf fc Henningl 
((1991 dust grains with thin ice mantles that coagulate 
for 1 5 yr at 10 6 cm~ 3 density (see Kau ffmann et all 
2008, Table A.l, for numerical values). This is con- 
sistent with previous studies of this subject, including 
sources as diverse as diffuse cl ouds and galaxies (see ref- 
erences in Bi anchi et al J 12003) . If one fixes the opacity 



use d for extinction observa t ions, i t thus appears that the 
the lOssenkopf fc Henningl (|1994l ) model opacities near 
1 mm wave length are too large by an average factor 
1.5«2.35 1 /2. 

In our emission-based column density estimates, 
we do therefore bas i cally adopt the aforementioned 
Osscnkopf & Hcnning (1994) opacities for dust emission 
observations near 1 mm wave length, but scale these 
down by a further factor 1.5 to bring mass estimates 
from extinction into harmony with those from dust emis- 
sion. As seen in Fig. [51 this procedure yields reasonable 
results, since the dust emission and extinction observa- 
tions for NGC1333 match within less than a factor 2 
(i.e., when comparing the masses of the most massive 
Bolocam-detected fragment to the one of the least mas- 
sive 2MASS-idcntified fragment). Incorrect assumptions 
about dust temperatures may cause most of this offset 
(see below). 

To some extent, correction factors between masses 
from dust emission and extinction are also influenced 
by spatial filtering affecting bolometer observations (Sec. 
14.2.41) . Future studies need to address this problem in 
more detail. Still, the aforementioned scaling factor 
aligns dust emission and extinction studies, which is the 
only aspect needed in our present study. 

4.2.3. Dust Temperatures 

For Perseus, IRosolowskv et al.l (|2008aD estimate gas 
temperatures between 9 K and 18 K. Assuming that 
dust and gas temperatures are similar, we do therefore 
adopt a temperature of (12.5 ± 2.5) K in our mass esti- 
mates. This temperature uncertainty results in a relative 
mass uncertainty of ~ 20% for Perseus. A few sources 
may have dust temperatures, and mass biases, outside of 
this range. 

4.2.4. Impact of Spatial Filtering 

Like all other upcoming, present, and past 
xound-based bolomete r-derived dust emission maps 
Kauffman n et al J 12008). the Bolocam maps of Perseus 
are not sensitive to structure larger than some 
instrument-dependent spatial scale. In its impact, 
this problem is similar to the spati al filtering in inter- 
ferometric imaging. For Bolocam, lEnoch et all (|2006l ) 
show the filtering scale to be of order V to 2' radius. 

Quantitatively, this removal of large-scale structure has 
an influence opposite to the impact of smoothing (Sec. 
14. ip : here, the relative mass-loss increases with spatial 
scale. This bias has to be considered carefully when 
using emission-based mass-size measurements for anal- 
ysis. Obviously, the true mass will be larger than the 
observed value. Filtering in bolometer maps is unfor- 
tunately too complex to be characte rized in a compact 
fashion; see [Kauffman n et all ((2008) for a few rules of 
thumb. For the particular case of the Bolocam maps of 
Perseus. lEnoch et all (|2006l) report losses < 10% for radii 
< 1', but do not sufficiently explore larger objects. 

4.3. Global Trends in Maximum Mass for given Size 

The mass-size tendencies seen, e.g., in Fig. [5] suggest 
to describe these trends with power laws of the form 

m(r) = mo (r / ra) b . (12) 
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Fig. 6. — Mass-size data for Perseus (solid lines). The data are from dust extinction maps probing large spatial scales, and Bolocam dust 
emission observations sensitive to small ones. The left panel highlights the nature of the data; das hed lines give the instrument-dependent 
sensitivity limits. Circles and the dotted line indicate calculation of the global mass- size slope (Sec. l4.3l . The middle panel highlights slopes 
derived by matching power-laws to parts of the data. Tangential slopes (Sec. 14.411 are derived in the same f ashio n, but for infinitesimal 
radius ranges. The right panel presents Perseus in the context of the reference mass-size relations from Sec. 13.11 Here, black solid lines 
highlight Bolocam-detected fragments in NGC1333, as well as the extinction-probed fragments containing this cluster. Other data are 
drawn using grey solid lines. 



In this, b is the slope of the relation, and mo is the in- 
tercept. As shown in Fig. [SI such laws trace, e.g., the 
mass-size relation of the most massive fragments in the 
0.04 < r/pc < 0.3 and 1 < r/pc < 8 radius intervals 
with maximum deviations of ±6%. 

Beyond such detailed descriptions of individual fea- 
tures, power-laws can be used to capture more global 
aspects of a cloud's structure. Consider, for example, the 
relation between the maximum fragment masses, m max , 
observed at radii of r sm = 0.05 pc and n g = 3.0 pc. 
(These radii are chosen to permit comparison with other 
clouds, as becomes more obvious below.) Based on these, 
we can define a global slope, 



frglob — 



ln[m max (rig ) /m max (r sm )] 
In [n g /r sm ] 



(13) 



As illustrated in Fig. [6l this slope is defined such that 
Eq. (|12[) connects the mass and size measurements for 
b = 6 g iob- In Perseus, 6 g i b = 1-21 ± 0.14, where we 
evaluate the uncertainty very conservatively by scaling 
the emission-based mass in Eq. (|13l) up and down by a 
factor 2. 

In Perseus, at the chosen scale of r = 0.05 pc (i.e., 
40"), the Bolocam dust emission map is not significantly 
affected by spatial filtering. Similarly, at r = 3 pc (40'), 
the extinction map is not biased by presence of stellar 
clusters. The measurements of 6 g i b are thus immune to 
these influences. 

4.4. Local Trends in Mass 

We now turn to slopes of infinitesimal tangents. As 
we explain below, the Bolocam data are not suited for 
this analysis, since they suffer from too strong spatial 
filtering. 

4.4.1. Method and Uncertainties 



In some cases, one may wish to use Eq. (|12|) to describe 
tangents to the data on spatial scales smaller than the 
one for which we calculate global slopes. This is demon- 
strated in Fig.|6]using tangents to the most massive frag- 
ments seen at given radius. In particular, it can be de- 
sirable to fit infinitesimal tangents to mass-size trends. 
For these, the slope reads 



b(r) 



d ln(m[r']) 



d ln(r' 



(14) 



Figure [T^a) shows slopes derived from mass and size dif- 
ferences between consecutive contours. As seen in the fig- 
ure, these data are rather noisy. We do therefore smooth 
the measurements. To do this, at radius r, we replace 
slope and radius by their respective unweighted arith- 
metic means, as derived within a smoothing kernel of 
width [r, 1.15 -r] (for a given cloud fragment, not permit- 
ting mergers). This yields the data shown in Fig. |7fb). 
Sometimes, however, the smoothing kernel is not filled 
well, and the data are still noisy. Thus, we do finally 
remove all data where the kernel is not filled to at least 
2/3. Figure [Tfc) shows this final result. 

The impact of noise can be estimated by propagating 
the mass uncertainties due to noise (Eq. ^ within the 
slope calculations. This yields 



^(ft) < 20 <r(m) 



(15) 



Because of the aforementioned smoothing, the numer- 
ical constant must be calibrate with our noise exper- 
iments (as done for Eq. [3]). To obtain an estimate 
of the expected uncertainties, we can repeat the mass 
and size substitutions done in the discussion of Eq. 
0. This gives maximum uncertainties ~ 60%. Since 
(r(b)/b oc (A^Ha) -1 r~ l (see discussion of Eq. [9]), the un- 
certainties are small for larger regions well above the de- 
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Fig. 7. — Calculation scheme for slopes (top to bottom) and im- 
pact of noise (black vs. yellow lines). We start with slopes directly 
calculated from mass and size differences for successive contours 
(panel a). These measurements are subsequently smoothed (panel 
b), and data are removed where the smoothing kernel is not filled 
(panel c). The black lines hold for the observed Perseus extinction 
data. As in Fig. [5] yellow and black lines indicate results from 
maps with and without artificial noise, respectively. The dotted 
line indicates the upper limit on slopes inherent to our method, 
6<2(Eq.0. 



tection threshold. In practice, uncertainties < 10% are a 
reasonable estimate for well-detected regions warranting 
detailed study. 

The slope difference due to smoothing is given by the 
first derivative of the mass bias due to smoothing (Eq. 
flO|) with respect to the radius. Including the usual cali- 
bration of numerical constants, we constrain the absolute 
smoothing-induced error to 



A6 < 0.4 to 0.6 



r k /r 



(1 + r k /r) 3 



(16) 



In a few cases, however, the error can be larger by a 
factor 2. In this paper, we reject regions with a diameter 
smaller than twice the beam diameter. Substitution of 
this limit (i.e., r^/r < 1/2) into Eq. (ITBl implies that 
slopes are overestimated by a number smaller than 0.1 
due to resolution. 

As mentioned in Sec. l4.2.4"l bolometer maps suffer from 
spatial filtering. This has an impact opposite to the in- 
fluence of smoothing, and artificially shallow slopes are 
measured from such maps. Since the filtering is very 
strong in Bolocam maps, we do not use these for the 
derivation of tangential slopes. 




1 
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Fig. 8. — Impact of limited spatial resolution on slope measure- 
ments. See Figs.[5{b) and [7] for explanations of mark-up. 



4.4.2. Results for Perseus 

Figures [3c) and |8] show tangential slopes for Perseus. 
The tangential slopes in the 1 < r/pc < 8 radius interval 
are in the range b — 1.2 to 1.7. At a given radius, the 
slopes for different fragments do often differ by more than 
their uncertainty. Also, in a given fragment the slope can 
change significantly with respect to radius. This means 
that it is not possible to describe an entire cloud by a sin- 
gle tangential slope. The observed tangential slopes bear 
no obvious relation to the global slope ^1.56 derived for 
the same radius range (Fig. [6]) . 

A slope of > 2 is actually not possible for a mass-size 
relation as defined by us, since this would require that 
the column density increases with radius (Eq. [3]). Since 
our search algorithm proceeds by decreasing the column 
density threshold, such fragments are not identified by 
our identification scheme. A slope b < would mean 
that the mass decreases with increasing radius. This is 
impossible in all but the most jolly insane models of cloud 
structure. 

5. SUMMARY & OUTLOOK 

5.1. Utility of and Outlook for Mass-Size Studies 

As we have shown above, our map characterization 
scheme yields reliable measurements of mass and size. 
From these, mass-size slopes and intercepts can be de- 
rived. Below, we describe how these data contribute to 
critical fields of star formation research. 

First of all, this approach permits a continuous char- 
acterization of cloud structure across a large range of 
spatial scales. This is just a desirable feature of any 
data analysis method, independent of the exact nature 
of the later analysis. The need for such a procedure 
led iRosolowskv et al.l (|2008bl ) to the development of the 
"dendrogram technique" . 

In star formation research, the basic mass-size mea- 
surements permit to compare fragment masses at a given 
spatial scale. Consider the classical case in order to see 
the advantage: usually, "cores" and "clumps" extracted 
from maps differ in their size. In this case, it is not clear 
what differences in masses mean, even if just a single 
cloud is considered. 

Spatially continuous cloud characterizations become 
particularly useful when comparing observations for dif- 
ferent molecular clouds. Usually, every cloud is stud- 
ied at a different physical resolution (i.e., pc per pixel). 
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In the classical case, mass measurements will thus usu- 
ally refer to vastly different spatial scales. With our 
approach, however, all scales are probed, and different 
clouds can be compared at the same physical scale. This 
is extensively employed in part II, where we study a sam- 
ple of clouds. 

The general utility of measurements of mass and size is 
known since long. For example, one can compare the ac- 
tual to virial masses or, more generally, masses predicted 
by theoretical cloud models. Equation (JHJ, for example, 
relates model fragment masses and gas temperatures. We 
shall not discuss such considerations here in detail. 

A property uniquely constrained by our method are 
mass-size slopes; these can only be measured via a scale- 
independent method. In particular, this gives access to 
the density structure of molecular clouds. For simple 
models of cloud structure, the mass-size slope is, e.g., 
directly related to the slope of the density profile (Eq.[7]). 
Such work on large-scale structure in molecular clouds is 
urgently needed, since cloud density profiles are presently 
not known on scales > 0.1 pc. Part II presents slope 
measurements for various molecular clouds, as well as 
several model mass-size relations. 

5.2. Summary 

This work studies the internal structure of molecular 
clouds by breaking individual cloud complexes up into 
several nested fragments. For these, we derive masses 
and sizes, as e.g. outlined in Fig. [T] Effectively, we per- 
form a "dendrogr am analysis" of a two-dim ensional map, 
as introduced bv iRosolowsky et"aTI (|2008bt ). 

The present paper establishes the method via a de- 
tailed analysis of the Perseus Molecular Cloud. Other 
solar neighborhood molecular clouds (< 500 pc; the Pipe 
Nebula, Taurus, Ophiuchus, and Orion) are discussed in 
the next paper of this series (part II). 

Power-laws of the form m(r) — mo (r/pc) b , with slope 
b and intercept mo, prove useful to quantify the rela- 
tions between mass, m, and size, r (i.e., the effective 
radius). Sections 14.31 and 14.41 discuss two different ap- 
proaches to define and measure the slope. We use global 



slopes to measure the relation between structure at small 
and large scales. This is done by connecting the mass-size 
measurements of the most massive fragments at small 
(0.05 pc) and large radius (3.0 pc) by a power-law (Eq. 
I13[) . Tangential slopes, on the other hand, are calculated 
infinitesimally at a given spatial scale (Eq. [T4| . The un- 
certainties in these properties are examined in Sees. 14.11 
and 14.21 (for mass and intercept), respectively Sees. 14.31 
and 14.41 (for slopes) . 

We conclude that our mass, slope, and intercept mea- 
surements provide a reliable method to characterize 
cloud structure. Our approach enables a continuous 
and reliable characterization of cloud structure in the 
0.05 < r/pc < 10 spatial range. This is not possi- 
ble using previous methods, since these are usually bi- 
ased towards a particular spatial scale (see, e.g., the 
CLUMPFIND analysis in Fig. 2]). Such comprehensive 
pictures of star-forming regions can be used to develop a 
more complete theoretical understanding of global cloud 
structure (Sec. 15. ip . 

A first observational and theoretical exploitation of 
this method is presented in part II of this series. We 
characterize, for example, the typical parameter space for 
solar neighborhood molecular clouds not forming mas- 
sive stars. Based on this, we chart a potential mass-size 
threshold for the formation of massive stars. Mass-size 
slopes are used to constrain large-scale density gradients 
within molecular clouds. 
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